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^\ ' A theoretical study is given of a new type of optical vortex in nonlinear anisotropic media. This 

is realized as a special solution of the two-component non-linear Schrodinger equation. The vortex 
►^ ' is inherent in the spin texture that is caused by an anisotropy of dielectric tensor, for which a 

Q I role of spin is played by the Stokes vector (or pseudo-spin). By using the effective Lagrangian for 

the pseudo-spin field, we give an explicit form for the vortex solution for the case of two types of 
optical anisotropy; that is, nonlinear counterpart of birefringence giving rise to the Faraday and 
Cotton- Mouton effects. We also examine the evolution equation of the new vortex with respect to 
the propagation direction. 
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Nonlinear optics has been a major subject in physics for long time |l|. The main interest is focussed on particular 
solutions of the basic non-linear equation that governs the light field or electromagnetic field. Among others, remark- 
able is an optical vortex, the existence of which has been early suggested in ||l| . Recently the detailed study has been 
(_i , carried out from both of theoretical and experimental point of view plpj . More recently experimental verification has 
2 ■ been also given for the multi- vortices H . The basic idea of the optical vortex follows an analogy with the superfluid 
vortex that is described by the complex order parameter for bose fluid [pj, namely, the equation for the light field, 
which is known as the nonlinear Schrodinger equation (NLS), is very similar to the Pitaevki equation. Thus it is 
2 , natural to expect the occurrence of the optical counterpart of the superfluid vortex. 
Q ' The purpose of this letter is to explore a possible new type of vortex in nonlinear and anisotropic media that are 

O , characterized by a variant of birefringence. We call this new type vortex an "optical spin vortex" . Our starting 
equation is the two-component NLS. The two-component NLS has been recently used for exploring an object of 
dynamic soliton H, which is the two modes induced waveguide leading to the soliton polarization dynamics. This 
^ ' work shares partly the basic idea with the present attempt in the point that two component NLS naturally incorporates 
^'' the polarization state of light. Indeed the concept of polarization plays an important role in modern optics, especially 
, _^ ' in crystal optics ||6|,[7|. The quantity describing the polarization state is realized by the Stokes parameters, which forms 
f^ , a pseudo-spin and is geometrically described by a point on the Poincare sphere. Thus the two component NLS can 
7-H ' be written in terms of the field of pseudo-spin, which is naturally achieved by introducing the effective "Lagrangian" 
0^ of fluid dynamical form. The use of the effective Lagrangian gives us a direct access to the analysis of the optical spin 
^^ ■ vortex. The main thrust is twofold: First we give an explicit form for the vortex solution by adopting some specific 
nonlinear birefringence; the concrete form is given by the nonlinear counterpart of the one causing the Faraday and 
Cotton-Mouton effects. This will be carried out by numerical way. We next examine the evolutional behavior of a 
vortex with respect to the direction of the wave propagation, 
i-rt ' Two component nonlinear Schrodinger equation. — First we derive the two-component NLS for the light wave 

^ ] traveling through anisotropic media. The procedure follows the one developed in a recent paper [p| . Suppose that the 
O • electromagnetic wave of wave vector k travels in the direction of z with the dielectric tensor e. The nonlinear nature 
of media implies that e has a field dependence in nonlinear form, the explicit form of which will be given later. Further 
we assume that e varies slowly compared with the wavenumber k. The z-axis is chosen as a principal axis of the 
dielectric tensor, namely, the axis corresponding to one of the eigenvalues of the dielectric tensor. In this geometry, 
; ^ ■ e is taken to be 2 x 2 matrix. Let us consider the field equation for the displacement field D, which is reduced from 
_C3 ■ the Maxwell equation |||,^: 

0.v^D.(^)%B.„ a, 

where V = I 4j, -4-] and {x, y) denotes the coordinate in the plane perpendicular to z axis. Now we put 

T>{x,y,z) = ^{x,y,z)e-x:p[iknQz] (2) 

with k = —, and no(= ^/co) means the refractive index for the case as if the medium is isotropic. The amplitude 
i{x, y, z) is written as f = *(/i, /2) — /lei -I- /2e2. We assume that f is slowly varying function of z besides {x, y), and 



ei and 62 denotes the basis of linear polarization. By substituting (||) into (|i|) and noting the slowly varying nature 
of f i.e., ^ mS^ fc|f|, we can derive the equation for the amplitude f, namely, we can only retain the first derivative 
II as well as the Laplacian with respect to {x,y) Q, hence 
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where A is the wavelength divided by 27r. This equation is regarded as a two-state Schrodinger equation where A just 
corresponds to the Planck constant and z plays a role of time variable. The components (/i, /2) couple each other to 
give rise to the change of polarization which is just the effect of birefringence governed by a 2 x 2 matrix "potential" 
1) = e — rig . V represents a deviation from the isotropic value and it becomes hermitian if the non-absorptive medium 
is concerned. From the hermiticity, the most general form of v is written as 



/3 — 17 vq — a 



(4) 



For later convenience, we transform the basis to the circular basis instead of the linear polarization (ei, 62), that is, 
e± = (l/\/2 ) (ei ± 162), which is written as (e+, e_) = T(ei,e2). Here T is given by 2 x 2 unitary matrix: 
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By introducing the wave function a,s tp ^ Ti = * (t/;^, "02)' ^^ have the Schrodinger equation for -0: 
with the transformed "Hamiltonian" 



(6) 



H = ThT-^ = V^ + V. 

no 
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The "field-dependent" potential V is written in terms of the Pauli spin; V — vq x 1 + X]i=i ^j'^'i- 

Effective Lagrangian for the pseudo-spin field. — We now introduce the "quantum" Lagrangian leading to the 
Schrodinger type equation, which is given by 



J = /"^t hx - H'j ipd'^xdz. 



(8) 



Indeed, the Dirac variation equation SI — recovers the Scrodinger equation. We write Lc — J ip''iX-4-ipd'^x and 
H = J ip^H'ipd'^x with {H' ^ T + V'), which are called the canonical term and the Hamiltonian term respectively. 
Here we note that V'{tp^ , V) differs from V in (7) and some relation holds between V and V', namely, V = V' + tp^ ^^ 
in order to recover the NLS. Having defined the Lagrangian for the two-component field ip, we rewrite this in terms 
of the Stokes parameters: This is defined as Si = ip^aiip, So = V'^IV' with i = x,y,z [p|jll|]. We see that the relation 
5o = 5^ + Sy + Si holds, namely, Sq gives the field strength; 5*0 = |D| . Using the spinor representation. 



V'l = VS'ocos-, ■02 



'So sin - exp[i(? 



(9) 



we have the polar form for the Stokes vector S — {Sx, Sy, S^) = {So sin 6 cos (f), 6*0 sin sin 0, So cos 9), which forms 
a pseudo-spin and is pictorially given by the point on the Poincare sphere. In terms of the angle variables, the 
Lagrangian is written as 



L = J^il^cos9)^^d'x-{Hr 
where the potential term V becomes 



V 



(10) 



F = 




J2visAdh 



(11) 



Here v'q and v'^'s are nonlinear functions of the field strength 5*0 as well as the angular functions {9, (p) and this feature 
may be required for a stability of special solution for the pseudo-spin field. The kinetic energy term Ht is given as a 
sum of three terms; Ht — ^ J VV'^VT/'d-^x = Hi + H where the first term becomes Hi ^ f j^^ (VS'q) d'^x, which 
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gives the energy that is needed for space modulation of the field strength, and the remaining terms are written as 



iJ = 



SqX 
m 



(V0)2+sin2-(V0)2[>d'x, 



(12) 



which is separated into two terms; H = H2 + H3, 

SqX- 



H2 = 
H3 = 



4no 
5*0 A 
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(13) 



Here if we define the "velocity field" v = (1 — cos0)V0, the first term is regarded as fluid kinetic energy inherent in 
spin structure, while the last term represents an intrinsic energy for the pseudo-spin which exactly coincides with a 
continuous Heisenberg spin chain |12[ . 

Vortex solution and its numerical evaluation. — We are now concerned with getting an explicit form for the specific 
type of solutions, namely, vortex solution for the two-component NLS. The solution we want here is a "static" solution, 
namely, we look for the solution that is independent of the variable z. For this purpose, we consider two types of 
anisotropy. 

(I) First we adopt the following nonlinear birefringence: 



gi^i-^i - i'2'4^2) 
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with the positive coupling constant g. This may be regarded as a nonlinear realization of birefringence that causes the 
Faraday effect. We have the potential V — J v'-^S^d^x = gSg J cos^ 9d^x. In what follows, we confine our argument 
to the case that Sq becomes constant. Physically, this corresponds to the constant background field with a proper 
core which is controlled by the profile of the angle functions {9, </>). A static solution for the one vortex is obtained 
by choosing the phase function <j> — ntan"-"^ i~)' with n = 1,2,- • •. being the winding number, together with the 
profile function 9 that is given as a function of the radial variable r(= ^/x^ + y^). Note that such a vortex becomes 
non-singular, namely, the velocity field v = (1 — cos 9) does not bear the singularity due to the behavior of 9{r) near 
the origin (see below). The static Hamiltonian is thus written in terms of the field 9{r): 
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where g' — _25o£2., xhe profile function 9{r) may be derived from the extremum of i?', namely, the Euler-Lagrange 
equation leads to 



(P9_ 

de 



ld9 



Ae 



■ sin 9 



■sin2e' = 



(16) 



where we adopt the scaling of the variable: ^ = \f^^- In order to examine the behavior of 0(^), we need a specific 
boundary condition at ^ = and ^ = 00. We impose 0(0) ~ 0, whereas at ^ = cx), there are two options: a) 9(po) = tt 
and b) 9{co) = 7r/2. If introducing the vector ni(^) = S/5o, we have 7713(0) = 1 for both cases a), b) and we have 
7713(00) = —1 for case a) and TO3(cxd) = for case b). This feature indicates that the pseudo-spin field which directs 
upward (left-handed circular polarization) at the origin changes to the state of dow7iwar(i( right handed polarization) or 
0M<war(i(linear polarization) with departing from the origin [see Fig.n^(a) and (b)]. We first consider the behavior near 
the origin ^ = 0, for which the differential equation behaves like the Bessel equation, so we see 9{S^) ~ Jn/2{0^ which 
satisfies 9{0) ~ 0. We examine the behavior at ^ = cx). This is simply performed by checking the stability for two 
cases mentioned above: (a) and (b). Now for the case (b), if putting 9{^) ~ ^ + a, with a the infinitesimal deviation, 
then we have the linearized equation a" — a ~ near ^ — 00, which results in a ~ exp[— ^]. This means that the 



solution with 9{oo) = f is stable. On the other hand, for the case (a) we have a" 



0, which gives a ~ exp[±i^] 



meaning the oscillatory behavior. This simply implies that the solution with 9(00) ~ tt does not converge to the 



stable solution which means that the case (a) is not relevant. Keeping mind of the above general feature, we here give 
a numerical solution of 6{S,) in Fig.||(b). 

(II) We next examine the nonlinear birefringence that is governed by the off-diagonal v matrix such that 

giipl->p2 



gir2^i ;■ (^^) 

This matrix may be regarded as a nonlinear counterpart of the birefringence that causes the so-called Cotton-Mouton 
effect 0. The corresponding potential energy becomes V = J X]i=i v[SiCpx = Sq J gicPx— Sq J gi cos^ 9(Px, the first 
of which is constant and should be discarded. Thus the resultant equation for the profile function 9{r) leads to 






2 ■ ,,, ^^,sin0--sin20^O (18) 



where the scaling variable S, = \f^\r with g'^ = _2i^lji. Que should note the "minus sign" in the last term. Due to 
this, if applying the same procedure in the case (I), we see that the behavior near the origin is given by ^(^) ~ -^n/2(C)i 
the modified Bessel function, which satisfies 0(0) — 0. At infinity, we get a stable solution for 6'(^) such that the 
boundary condition 0(oo) = tt is satisfied. This feature is opposite to the previous case, namely, the solution satisfying 
the boundary condition 0(oo) = ^ oscillates so it should be omitted. The numerical result is also given in Fig||(a). 

In summary, the vortex solutions in these two cases show up quite different behaviors each other, due to the 
difference of nonlinear birefringence. From Fig. 2 (a) and (b), we can estimate the scaled vortex-core size ^c from 
^c0(oo) — S where the area S surrounded by the solution curve and the asymptotic line 9 = 6{oo) — n (for case 
(a)) or 7r/2 (for case (b)). The ^c is just the mean value for the area S. Using the numerical results, we get 
^c = S/^ = 0.714 for (b) and ^c = S/tt = 0.61 for (a), which are consistent with half the scaled coherent (healing) 
length ^cohf^ = rcoh/Vg' = l/\/2 ~ 0.7 The real (unsealed) core-size Vc is given by Tc = ^c/Vg' for both (a) and (b). 
If the characteristic wavelength of light is larger than the core size; the core should be detectable: its condition is 
given by ^^ < A^F- 

Evolution equation for vortex. — Having demonstrated the explicit form for the vortex solution, we now consider 
the evolutional behavior for a single vortex with respect to the propagation direction z. Following the procedure used 
in the magnetic vortex |1J], let us introduce the coordinate of the center of vortex, R(2:) = {X{z), Y{z)), by which the 
vortex solution is parameterized such that 0(x — 'R{z)) and 0(x— R(z)). By using this parametrization, the canonical 
term Lc, the first term in ([l(]|), is written as 

Lc = ^ fvTld^x (19) 



where we have used the relation: ^ = ^R, 4^ = ~V4> with R = ^. This can be obtained by the "Euler-Lagrange" 
equation for R, which gives the "balance of forces" 



d dLc dLc dH 



Using eq.(|l9|), we get 
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where the k is the unit vector perpendicular to the xy-plane. Here a is defined as 

In deriving (pl|), we have used the relation ^^ = ~^^- '^^'^ integrand of a is nothing but the vorticity which 
we put oj. Using the expression for the velocity field in Eq.dlfl), we can write ou in terms of the angular functions: 
w = (V X v)^ = sin6'(V0 x V0), or in terms of the spin field m 

/„ N /dm dni\ ,„„, 

(Vxv)..m.(^— x_j. (23) 



The equation (E3f) is an optical counterpart of a topological invariant of hydrodynamical origin |14| , which is written 
as 

cr = /sin 0de A d(l) (24) 

Js 

where S stands for the area in the pseudo-spin space {9,(f>). a has a topological meaning, which depends on the 
boundary condition for 0{r). Namely, for the case a) corresponding to the boundary condition 6{oo) = tt, the vortex 
configuration gives the mapping from the compactified two-dimensional space i?^ U oo ~ S*^ to the pseudo-spin space 
82- Hence a in d24) has a meaning of the degree of mapping for S2 -^ S2 leading to the topological invariant 
a: a = n (n^integer). For the case b) corresponding to the boundary condition 6(r) = ^, the mapping becomes 
^2 — > S'2/2(hemisphere), so we have the topological invariant a — njl. The appearance of such two types of topological 
invariant is characteristics of a new type vortex presented here. 
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FIG. 1. The profile of the non-singular vortex; a) The 
case of 6{oo) — n and b) The case of ^(oo) = |-. 



FIG. 2. The profile of the function 9{^) for two cases of 
boundary conditions: (a) 9{oo) = n and (b) 9{oo) = ^. 



